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Abstract 

A general class of solutions of Einstein's equation for a slowly rotating 
fluid source, with supporting internal pressure, is matched using Lich- 
nerowicz junction conditions, to the Kerr metric up to and including first 
order terms in angular speed parameter. It is shown that the match ap- 
plies to any previously known non-rotating fluid source made to rotate 
slowly for which a zero pressure boundary surface exists. The method 
is applied to the dust source of Robertson- Walker and in outline to an 
interior solution due to McVittie describing gravitational collapse. The 
applicability of the method to additional examples is transparent. The 
differential angular velocity of the rotating systems is determined and the 
induced rotation of local inertial frame is exhibited. 



1 Introduction 

In the period since the discovery of the Kerr [1] metric which describes ana- 
lytically, the asymptotically flat, vacuum gravitational field outside a rotating 
source in terms of Einstein's field equations, there have been many attempts 
to find closed interior solutions which match the exterior smoothly. In general 
terms attempts to find solutions have proved unsuccessful as has been described 
by Pichon and Lynden-Bell [2]. One difficulty has been the considerable math- 
ematical complexity in solving Einstein's equations, see for example, Krasinski 
[3], Chinea & Gonzalez- Romero [4]. This has led to an 'embarrassing hiatus' 
according to Bradley et al [5] in the number of potential interior solutions avail- 
able for matching which in turn has contributed to a lack in the development in 
the theory of differentially rotating fiuid bodies in general relativity. Even for 
the case of the remarkable and much quoted, Wahlquist [6] closed form interior 

'Present address: Department of Physics & Astronomy, Cardiff University, Queen's Build- 
ings, 5 The Parade, PO Box 913, Cardiff Cf24 3YB. email: Ron.Wiltshire@astro.cf.ac.uk. 



1 



there is no possible fit to the Kerr exterior as has been recently been shown by 
Bradley et al [5]. Only for the important case of thin super-massive rotating 
discs, supported by internal pressure have analytic sources for the Kerr metric 
been found (Pichon and Lynden-Bell [2] ) . Yet it is important to develop further 
the relativistic theory of rotation since it has considerable potential application 
in astrophysics, for example, in the description of the gravitational collapse of 
rotating matter, quasars, or potential sources for gravitational radiation. 

To this end one way forward is perhaps to adopt a linearised perturbation 
technique wherein, in the first instance, the interior source is rotating only very 
slowly. Following the successful match to the Kerr exterior one would proceed 
to develop higher order perturbation methods to describe sources with higher 
angular velocity. The use of perturbation techniques in general relativity are of 
course common. In the context of rotation they have been applied successfully 
by Hartle [7] to equilibrium configurations of cold stars. In the case of non- 
equilibrium configurations Kegeles [8] has applied the method to Robertson- 
Walker dust sources up to the first order in angular velocity parameter although, 
the results are somewhat restrictive and are not suitable for application to 
sources supported by internal pressure. 

The apppropriate junction conditions for solutions of Einstein's equations 
are extensively considered in the literature for example, Misner et al [9], Mars 
& Senovilla [10], Stephani [11], Hernandez-Pastora et al [12] where the main 
focus of attention concerns the methods of Darmois [13] and Lichnerowicz [14] 
which have been shown to be equivalent by Bonnor and Vickers [15]. In the 
Darmois approach it is necessary that the components of the metric tensor, 
and also the extrinsic curvature for the Kerr exterior and the interior source 
are continuous at the boundary surface. A common coordinate description 
of source and exterior is not required. On the other hand, the Lichnerowicz 
approach requires a common 'admissible' coordinate system wherein both the 
metric tensor and its first partial derivatives are continuous at the boundary 
surface. Although flexibility and covariance has ensured the extensive use of the 
Darmois approach in the literature, the requirement to define the 'admissible' 
system by Lichnerowicz has proven to be of benefit for the analysis below. 

The main aim here is therefore to develop the perturbation method for 
rotating systems up to and including first order terms in the angular velocity 
parameter and to find general solutions of Einstein's equations for perfect fluid 
bodies which fit the Kerr solution. The solution of Einstein's equations for 
the non-rotating system will be assumed given as for example, the comoving 
cases presented by McVittie [15], Kustaanheimo [16], Bonnor & Faulkes [17], 
Chakravarty ct al [18] and other described by Kramer et al [19]. 

Thus in the following the spherically symmetric source will be a given solu- 
tion of Einstein's equations described by means of the metric 

dn'^ = dO"^ + sin^ ed<p'^ , (1) 

where A = A(^, ry), fi = /v, (^,77), r = r{^,ri) .In the following the components 
of this metric will denoted by the tensor gab- The source boundary surface will 
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be described by the equation: 



ib = Fiv) , (2) 

where F (rj) is some function of rj alone and where the suffix '6' will denote 
evaluation at the boundary throughout. Lichnerowicz conditions will be used 
to join the interior to the Schwarzchild exterior which will be written in the 
form: 

dn^ = S^ + sin'^ed^^ , (3) 

where: 

e^ = l-^>0. (4) 

It will be assumed that the coordinate description 9, (j), rj) is a suitable ad- 
missible system where the boundary conditions apply and that this is related 
to the (S, ^, ^, n) description of the exterior by means of the transformation: 

S = S(C,77) , ~9 = 9 , ~^ = 4> , n = n(e,ry) . (5) 

It follows that (^) can be transformed to the form: 

da^ = (e^n^ - e-^S^) cir?2 + 2 (e^n^n^ - e-^S^S^) dCdjy 

- (e-^S| - e^n|) df - S2 (e, r/) dn^ (6) 

where the suffices ^ and rj mean partial derivatives with respect to ^ and r/ 
respectively. The metric components of (|6|) will be denoted by 7^^. In addition 
units are chosen such that c = 1 = G. 

The approach adopted here therefore will be at first to apply the Lichnerow- 
icz junction conditions, in which the components of Qab and, their first partial 
derivatives describing a spherically symmetric non stationary fluid sphere is 
matched continuously to and, their first partial derivatives describing the 
Schwarzchild metric. The conditions will then be considered in the context of 
slowly rotating systems. The results are used to construct new solutions of 
Einstein's equations which are applicable for slow rotation and as an example 
a rotating dust source is considered and the results of Kegeles [8]. In addition 
in a further example the McVittie [15] source is 'set' into slow rotation. 

2 Application of the Lichnerowicz junction condi- 
tions 

By means of (||) and (^) the continuity of the metric components g^i and 7^,; , 
kl = 22 and first partial derivatives across the boundary imply that 

Sfe = rft , {SJ^ = {rj^ , {S^}^ = {r^}^ (7) 
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where, for any function X = X {£^,rj): 

I OK J i^=F{ri) I or] J ^^=F(ri) 



(8) 



However, for any X^, : 

dXb 
drj 



d{X^}, 



b 



{Xrinlb - +^vi^(Ub 

and so applying this to {E^^}^ and it also follows that: 

{^v^}b = {^'(Jb + {'^^Sife - ^'7 {%}& ' 

{S,,}, = {r,,},-F2{r^J, + F2{E5j, . (10) 



Thus equations (|10D may be used to determine {S^^}^ and {S^^}^ once F and 
{S^^}^ have been determined. In particular, when = then = 
and {S^^ljj = {rriri}fy- 

Furthermore from (||) and @ continuity of the component Qki and 7^^, 
fe/ = 44, 14, 11 of the metric across the boundary gives rise to three further 
relationships that: 



{n.}, = {e-^(e^^ + e-^E^)} 



b 



mb = {e-'^ie-'^n-e'^m . (11) 



and the physical restriction: 

|e2A5.2 _ g2(M+A)+^ _ g2M5^2|^ ^ Q ^ (12) 

where: 

Equation ( |T2| ) is merely the condition: 

{m(^,r?)}^ = m . (14) 

where the mass function m (^, r]) is as usual defined, in terms of the Riemann 
tensor through: 

^ (e, r?) = ^ijL = ^ (l + e-^V^ - e-^/^rf) . (15) 
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This condition may be used to simplify (|TT|) with the result that: 



{n.}, = {e-^e^-'^s^}^ , {na, = {. 



(16) 



Consider now the continuity of the partial derivatives of the metric com- 
ponents Qki, and 7;.^ for kl = 44, 11 across the boundary surface. It follows 
immediately that: 



d_ 
drj 
d_ 



e"- I e^^ + e 



(17) 



and: 



d_ 



J J fe 



' 6 

(18) 



Furthermore direct expansion of the consistency relation {911^/9^}^^ = {dll^/drj}j^ 
in (ITSI) gives rise to: 



{Er,^- ^i^r^- X^r^}^ = 



(19) 



and, using (10) may be expanded in terms of the Einstein tensor component 
Gf, to give: 



+ {r^^ - E^a, = . 



(20) 



Furthermore using (0) , ( [l8|) and it follows that the continuity of the partial 
derivatives of the metric components g^i, and for kl = 14 can also be written 
as: 







d_ 

dx 



(21) 



where x = i], ^. With the aid of (|To|), these conditions may also be written in 
terms of the Einstein tensor components Gf and Gf . Thus: 



F^{Gtr^-Giri:}^ + {Glr^-G\r,}^ = . 



(22) 



The second equation in ( p2D may also be expressed in terms of the mass function 
to give: 



dm ^ 
drj ji 







(23) 
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The equations (pO|), and ( |2^ ) may be substituted one into the other to obtain 
the following simplification: 

{Gi-F,G^}, = , {G1-F,GO, = (24) 

{re2^Gl}, + 2{S55-r^a, = 0. (25) 

These equations may be used to determine F, {S^^}^ and also express the 
restriction (p3|). 

Note that the above results have been obtained without reference to the 
nature of the energy momentum tensor . However if it is now supposed that 
the source is a perfect fluid then (p^, (p3|) become: 



Pb = , {u^-F^u^}^ = (26) 

as expected, where p t]) is the source pressure and are the components of 
the velocity four-vector. Further in a comoving description where = then 
F {ij) is a constant, again as expected. 

Thus direct application of the boundary conditions specifies conditions or 
limitations on the functions T, (,^, rj) and 11 (^, r]) but in no way defines them 
uniquely. However, the nature of these transformation functions must be estab- 
lished to complete the definition of the admissible system. It is straightforward 
to check that the transformation functions: 



ra=3 



(27) 



(28) 



S(^,r/) =r(^,r/) 

and 

n (e, r?) = n, + {Hg}, (C - F) + ^"^^^^ (g - FT 

n=3 

with 

U, = J ({n,}, + F,{njJdr? , (29) 

and Dn = Dn {rj), En = En{r]), n > 3, are consistent with equations (0), (p^), 
(|l|), (0), (|l|) and (H). The ambiguity of S (^, r?) and H (^, r/) is clear through 
the arbritrary definition of Dn and En- 



3 Junction conditions for slowly rotating systems 

Consider now a general first order rotating source given by 

dcT^ = e^^dr/^ - e'^'^df - r^dn'^ - 2r^ sin^ {9) q {Yd^dcp + Xd^drj) (30) 

where X = X(^,r/), Y = 1^(^,7?) and g is a small angular speed parameter 
whose square terms and higher are negligible. Note that up to and including 
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this order the source boundary may still be written as ^ = .^f,. The exterior 
metric now the Kerr solution which is transformed to the form: 



da^ = (e^n^ - e-^S^) + 2 (e^n^n^ - e-^S^S^) d^djy 
- (e-^S| - e^n|) de - S2 (e, rj) dQ^ 

— n^dC^ + —Ur^d(f)dr] j (31) 

where, a is the angular speed parameter with negligibly small quadratic terms.. 
The additional continuity conditions applied to X = X (^, rj) and Y = Y {(^, r/) 
are obtained by comparing (30) with (31) so that with a = q: 

Moreover from (^) the continuity conditions applied to the derivative of 
the metric tensor require: 

W.^{|(^)}, . (33) 

and: 

and so using (|l5|): 

(35) 



4 Solution of Einstein's equations for slowly rotating 
systems 

It will now be shown that solutions of Einstein's equations satisfying the junc- 
tion conditions do exist. Suppose that Einstein's equations for a perfect fluid 
source are written in the form: 

Gt = -87TT,- , T^ = (p + p)u-u,-5tp , (36) 

where p, p are the source density and pressure and u"" are the components of 
the velocity four-vector with the property that u°'Ua = 1 and further suppose, 
for simplicity that comoving fluid spheres satisfying G\ = and G\ = are 
given for a non-rotating source. In addition, it is noted that the source density 
and supporting internal pressure for the metric ( ^0|) are not affected by the 
addition of a rotation speed parameter up to and including order q. However, 
for the slowly rotating systems it is necessary to ensure that the components of 
the velocity four- vector satisfy = u^. Now it is straightforward to show that 
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Gl = = Gl ,Gl = = Gl identically for the metric (|3^), whilst Gl = = Gf 
will be satisfied so long as: 

{Y,,-X^,) = {Y,-X^)(^X, + t,^-^^ . (37) 

Thus Einstein's equations are satisfied whenever: 

y, = x^ + ^^%^ (38) 



where h is an arbitrary function of ^. So comparing ( |35D with (|38| ) it is seen 
that the boundary conditions are consistent with Einstein's equations provided 
that: 

hb = 6fh. (39) 
In addition, the angular velocity L (^, rj) of the source is given by: 

and, a particle moving in the field of (^) will have zero angular momentum 
whenever M3 = 0, so that the quantity: 

ng ^ gsin^ ge~^-^/tg ^^^^ 
U4 IQirr'^ {p + p) 

will also be zero for such a particle. It follows that the induced angular velocity 
0, (^, ry) of the inertial frame is given by: 

n{^,v) = -qX , (42) 

and that with (|l|) and (|3|): 

This confirms that the 'frame dragging' effect decreases inversely with the cube 
of r as is well known (for example, Schutz [20]). 

Clearly it is now necessary to determine those functions X (^, ry) and Y (^, rj) 
satisfying the junction conditions which also satisfy (|38[). These are established 
by writing: 

XiC,v) = Xfe + (e-e,){Xj, + ^(^,r7) , 
Y{tv) = n + ^(e,r/) , (44) 
where {^I'5}J = and using (^) and ( p9|) in (H) gives rise to: 



$(e,r/) 



+ ^^^^^ J >> 



dr] . (45) 



Thus ( ^4|) with (^) is solution for X (^, r;) in terms of Y (^, rj) and satisfying 
the necessary boundary conditions (|33|) to (|35|). 
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5 A slowly rotating dust cloud 



Suppose that a slowly rotating dust cloud is described by a perturbed Robertson 
Walker metric of the type: 



da' 



dri' - R' (r/) 



de 



1 - 



+ i'de' + r sin^ {0) d^' 



2^' sin^ (9) R'q [Yd^d^ + Xd^dr]) 



(46) 



It is easy to see that {m ry)}^ = m from ( p^ ) and pb = from ( p6[ ) give rise 
to: 



2m 



k 



The boundary relations for S (^,?7) are: 



(47) 



m 



3m 



ebii-kii) 

whilst for n (^, Vj) : 



2m 



(48) 



I drj 

I- 



{n.}, = {e-^},(i-fee 



e-^(l- A:^')^^^ 



1^ 




1 


1 drj 


.(1-fce')^ 





N I v2 

e I 2.^ 



(49) 



Furthermore the final solution (^ ) with (^ ) may be written as: 

d ( 2me-^(l - kf)hT.^^ 



= — — 4^^T^ — ^ + (e-e 



with: 



+*(e,r?) , 

2mfi^{e-^}^ 

an-kiD'^R'' 

^(e,r/)= h^id'n + 



MO 



6m 



(50) 



(51) 



This is a general representation of a collapsing, slowly rotating, Robertson- 
Walker dust cloud and generalises the results given by Kegeles (1978) who does 
not give transparent forms X rf) and Y (^, rj) . 
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6 Example with non-zero internal pressure 



Consider now, in outline a representative of a broad class of solutions in which 
the rotating source has non zero supporting pressure. The chosen source is the 
McVittie [15] solution although any spherically symmetric, comoving solution 
of Einstein's equations with an added rotation term could have been used. The 
only restriction on their use is that any solution must have a boundary surface 
where the mass function is constant, or, equivalently the supporting pressure is 
zero. 



In this case consider (pO) with: 







fm 






fm 


p 


2^S 



e 



1 _ 

1 + 



, fiO = {l-ciep (52) 

where S = S (rj) and ci is constant. The respective internal density and sup- 
porting pressure are given by:: 

K^S^ + 128^3^5 ff2 
8-P = J^s^ ■ (53) 



2K^SS„„ + U^K^SSl - bK^Sl + 256^^5^ if - I) 
SVTD = — f54) 

where K = 2^S+rhf . Since this metric is described by a comoving observer then 
-F,^ = and the boundary surface is given by the constant ^ = ^j,.This surface 
is defined through the constancy of the mass function ( |l^ at the boundary or 
by differentiating this with respect to r] to obtain the zero pressure boundary 
condition (|2^). Using the equations one may calculate: 

^^^\^^2[esH2iS-f-K)]i\^ ^^^^ 

The transformation functions, Ef,, {Sg}^, {S^,{Sg^}j^, {E^^}^ and {S^^}^ 
may be calculated directly using (0), (|lO|) and (|^) whilst , lib, {n^};,) {nr?};,, 
{n^jj}^, {I[jjn}fy and {11^^}^ may be determined easily from ([l^), (^) and (p^). 
Some simplification of the results will occur through the use of (|l5|), (|55[). The 
resulting somewhat lengthy expressions need not be reproduced here. 

To obtain the corresponding rotating solution it is necessary to calculate 
each of X^,, {X^}^ ,{X^}j,yfe,{yg}^ and {i^r?};, to ensure a smooth match to 
empty space-time of which Xb ,{X^}f^ ,Yb, and {Y^} are also required for the 
explicit determination of the solution (^) and (|45|). These calculations are 
straightforward. 
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Also note that to obtain the final solution it is necessary to evaluate the 
integral (|45| ) for which the integrand contain terms of the type: 



— 7^ — = ~. — } TT^ ^^^^ 

es^ (i + 

where , /15 = 6m. Clearly the resulting solution will therefore not have a closed 
form. None the less the solution for a rotating version of the McVittie (1933) 
has been established formally, although it is interesting to note that a solution 
so elegant in the absence of rotation has such an awkward form when slow 
rotation is included. 



7 Conclusion 

Solutions of Einstein's equations for slowly rotating time varying sources sup- 
ported by internal pressure have been presented. It has been shown that any 
known collapsing, or expanding, fluid source known by a comoving observer 
not to be rotating may be 'made' to rotate slowly and also matched smoothly 
to the Kerr exterior at all times provided that a zero pressure boundary sur- 
face exists. In each case the source rotates with an angular velocity which is 
inversely proportional to the sum of its internal density and pressure and the 
induced rotation of the inertial frame of reference is explicitly clear. 

It has been shown that rotating solution may be expressed in terms of 
the expression J e'^~^^^r~'^dr] and thus it follows that this must be integrable in 
closed form for the existence of analytic solutions. It is interesting to speculate 
that solutions of Einstein's equations have previously been found without the 
need to give the nature of this expression any consideration. For example 
it has previously been common to express comoving systems in terms of the 
isotropic coordinate r = ^e^ which may not necessarily be an appropriate choice 
for the development of rotating descriptions. This matter is currently under 
investigation. 

The analysis here has been presented here with accuracy up to to an in- 
cluding first order terms in the angular speed parameter. The natural extension 
of the work to the matching of second order case with Kerr space-time, where 
the source boundary is more complex in nature, will be reported in the near 
future. Extensions of this approach to higher orders may result in an analytic 
description of sources supported by internal pressure and which are surrounded 
by empty space-time. 

Finally, although the focus here has been upon slowly rotating compact 
bodies it is interesting to note that the Robertson- Walker metric with the added 



rotation term (46) may be considered as a cosmological model. In this case the 
junctions conditions need not be applied. The expressions for pressure and 
density remain homogeneous and the angular velocity term ( ^0|) may also be 
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made independent of ^ by choosing: 




(57) 



{I -key 



and X = X (rj) . The solution of Einstein's equations (^) then may be written 
as: 
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